Witnessing nonclassical multipartite states 
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We investigate a witness for nonclassical multipartite states based on their disturbance under local 
measurements. The witness operator provides a sufficient condition for nonclassicality that coincides 
with a nonvanishing global quantum discord, but it does not demand an extremization procedure. 
Moreover, for the case of ^-symmetric systems, we rewrite the witness in terms of correlation 
functions so that classicality is found to necessarily require either vanishing magnetization in the 
invariant axis or isotropy of the two-point function in the transverse spin plane. We illustrate 
our results in quantum spin chains, where a characterization of factorized ground states (with 
spontaneously broken Z2 symmetry) is achieved. As a by-product, the witness will also be shown to 
indicate a second-order quantum phase transitions, which will be illustrated both for the XY and 
Ashkin- Teller spin chains. 

PACS numbers: 03.65.Ud, 03.67.Mn, 75.10.Jm 



I. INTRODUCTION 

The characterization of quantum correlated multipar- 
tite systems has been a central topic in quantum infor- 
mation theory (see, e.g., Refs. |l|-|3[ for recent discus- 
sions). Such systems are described by quantum states 
that cannot be represented by joint probability distribu- 
tions. Remarkably, the existence of quantum correlations 
may occur even in the absence of entanglement and still 
be a resource for a number of quantum information ap- 
plications (see, e.g., Refs. 0-[l])- Measures of quantum 
correlations provide then a valuable tool to recognize 
useful quantum states for a certain task. In this con- 
text, an important measure is given by quantum discord 
(QD) [9(. It arises as a difference between two expres- 
sions for the total correlation in a bipartite system (as 
measured by the mutual information), which are classi- 
cally equivalent but distinct in the quantum regime. QD 
has been generalized to the multipartite scenario by dif- 
ferent approaches [Iol-[l3j|. Unfortunately, both bipartite 
and multipartite versions of QD involve extremization 
procedures, which makes their computation viable only 
for very small systems or for particular classes of states. 

In this context, a possible strategy to identify com- 
posite systems that do not exhibit purely classical cor- 
relations is to look for a witness operator, namely, an 
observable able to detect the presence of nonclassical 
states. A witness W is denned here as a Hermitian op- 
erator whose norm ||W|| is such that, for all nonclassi- 
cal states, ||W|| > 0, with a convenient norm measure 
adopted. Therefore, ||W|| > is a sufficient condition 
for nonclassicality (or, equivalently, a necessary condi- 
tion for classicality). Several investigations for estab- 
lishing conditions to witness bipartite nonclassical states 
have recently been carried out [T^ - [l7j . Our aim here 
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is to provide a systematic approach to obtaining a wit- 
ness for nonclassicality in multipartite states, which will 
be based on their disturbance under local measurements. 
In turn, such a witness will naturally detect states with 
a nonvanishing multipartite QD. As an illustration, we 
will consider arbitrary .^-symmetric states, where the 
witness operator can be conveniently rewritten in terms 
of correlation functions so that classicality is found to re- 
quire either vanishing magnetization in the Z direction or 
isotropy of the two-point function in the XY spin plane. 
In particular, we will apply such results in quantum spin 
chains, where it will be shown that ||W|| exhibits a com- 
pletely different behavior in the cases of ^-symmetric 
states and states with spontaneous symmetry breaking 
(SSB). In the latter characterization of nontrivial 

factorized ground states (with Z2 SSB) will be achieved. 
As a by-product, the witness will also be shown to ex- 
hibit a nonanalytical behavior at second-order quantum 
phase transitions. This nonanalyticity will be revealed 
by the first derivative of ||W|| and will be used as a tool 
to characterize second-order quantum critical points in 
both the XY and Ashkin- Teller spin chains. 



II. WITNESS FOR NONCLASSICAL STATES 

Let us begin by considering a composite system in an 
^-partite Hilbert space rl = %Ai ® Ha 2 ® ■■■ ® Wa n ■ 
The system is characterized by quantum states described 
by density operators p € B(rl), where B(rl) is the set 
of bound, positive-semidefinite operators acting on rl 
with trace given by Tr p = 1 . In order to define a 
multipartite classical state, we will generalize the idea 
of nondisturbance under projective measurements intro- 
duced for bipartite states in Ref. [HJ]. Indeed, we will 
denote a set of local von Neumann measurements as 
{IT.} = {11^®- ■ ■ U 1 / }, with j denoting the index string 
(*!"■• *at)- Then, after a non-selective measurement, the 
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density operator p becomes 



(i) 



This operation can then be used to define a classical state. 

Definition. If there exists any measurement {IX,-} such 
that $(/o) = p then p describes a classical state under 
von Neumann local measurements. 

Therefore, it is always possible to find out a local mea- 
surement basis such that a classical state p is kept undis- 
turbed. In this case, we will denote p <E C N , where C N is 
the set of A-partite classical states. Observe that such 
a definition of classicality coincides with the vanishing of 
global QD Note also that $(p) by itself is a classical 
state for any p, namely, ^(^(p)) = 3>(p). Hence, <&(p) 
can be interpreted as a decohered version of p induced 
by measurement. A witness for nonclassical states 

can be directly obtained from the observation that the 
elements of the set {TEj } are eigenprojectors of p. This 
can be shown by the theorem below (see also Ref. [l8| V 



Theorem 1 p € C 



N 



[p,Uj] = (Vj), with Tlj = 



rP 1 



IT/ 

A N 



and j denoting the index string (ii ■ ■ ■ in)- 



Proof. If p € C N then = p. Then, similarly as in 
Ref. [l8[, a direct evaluation of J^j [P^j] [p^j] yields 
J2j [p,~Rj] [p,^ = $(P 2 ) - P 2 - However, $(p) = p also 
implies that $(/9 2 ) = p 2 . Therefore, 

£>,iy [ P ,u^ = o. (2) 



Observe that, given a composite multipartite state 
p, it is rather simple to evaluate the commutator 
[p, pa 1 ® • • • ® Pa n ], with no extremization procedure as 
usually required by QD computation. From the neces- 
sary condition for classical states above, we can define a 
witness for nonclassicality by the norm of the operator 
[p, pa,--- PA N }- Indeed, if p € C N =*> ||W|| = 0, where 



W= \p,p Al 



® PA„ 



(5) 



Therefore, ||W|| > is sufficient for nonclassicality. 

For concreteness, we will take ||W| as defined by the 
trace norm, namely, 



|| W)| = TrVWWt 



(6) 



The witness ||W|| provides a sufficient - but not neces- 
sary - condition for nonclassicality. As a simple example, 
consider a two-qubit system AB prepared in a Bell state 
\tp) = (jOO) + |11))/V2, where {|0), |1)} denotes the com- 
putational basis for both qubits A and B. The reduced 
density operators for each qubit are given by maximally 
mixed states, with pa <8> Pb proportional to the identity 
operator. Therefore, by taking p = we have that 

[p, pA® Pb] = 0, even though p describes a maximally 
entangled state. However, as we will show below, such 
a witness can be a very useful tool to investigate non- 
classicality in other (non-maximally entangled) classes of 
quantum correlated states in many-body systems. 



Hence, p G C N =$> [p,Tlj] = 0. On the other hand, if 
[p, Uj] = then {n,-} provides a basis of eigenprojectors 
of p. Then, from the spectral decomposition, we obtain 
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Hence 



which immediately implies that &(p) = 
[p,Uj] =0=^ P eC N . u 

We can now propose a necessary condition to be 
obeyed for arbitrary multipartite classical states. 

Theorem 2 Let p be a classical state and pAi the re- 
duced density operator for the subsystem A{. Then 
[p,PA! ®---®Pa n ] =0. 

Proof. From theorem [T] if p G C N then the spectral 
decomposition of p yields Eq. ([3]). Therefore, 



4 = E^ n x- 



Hence, 

[P, PA 1 C 



PA 



by direct evaluation, 
•••(g) pa n ] = 0. ■ 



(4) 

we obtain that 



III. NONCLASSICAL STATES IN BIPARTITE 
SYSTEMS WITH Z 2 SYMMETRY 



As an application of the nonclassicality witness, let us 
translate the witness operator ||W|| in terms of correla- 
tion functions and magnetization in quantum spin sys- 
tems. To this aim, we will consider an interacting pair 
of spins- 1/2 particles in a chain, which is governed by an 
arbitrary Hamiltonian H that is both real and exhibits 
Z2 symmetry, i.e. invariance under 7r-rotation around a 
given spin axis. By taking this spin axis as the z direc- 
tion, this implies the commutation of H with the parity 
operator 0^ 1 <rf, where N denotes the total number 
of spins and of is the Pauli operator along the z-axis 
at site i. From this symmetry, the two-spin reduced 
density matrix at sites labelled by i and j in the basis 
{ I ft), I t0, 1 It), I 44)}, with I t) and | |) denoting the 
eigenstates of c 3 , will be given by an A-state, reading 



Pab = 



( a / 

61 z 

z b 2 

V/ d 



(7) 
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In terms of spin correlation functions, these elements can 
be written as 

a = \{ l + G l z + G{ + G%) , z = \ (G« + G« ) , 
h = \ (1 + G\ - G{ - G% ) , / l(G£ G%) , 
h = \ (1 - Gl + Gi - G« ) , 

d = ~ (1 - Gl - G{ + G»,) , 



(8) 



where G^ = (c 2 ) (fc = i, j) is the magnetization density 
at site k and G^ = (c^ct^) (a,/3 = x,y,z) denote two- 
point spin-spin functions at sites i and j, with the expec- 
tation value taken over the quantum state of the system. 
Note that, in case of translation invariance, which will 
be assumed here for simplicity, G Z =G Z = G 2 (V fc, k') 
and, therefore, b\ = &2 = b. In order to evaluate the 
witness in terms of the magnetization density and spin- 
spin correlation functions, we compute the commutator 
W = [p, pa <8> pb], yielding 



w = 



k 


-k 



(9) 



with k = f (b + df - (a + b) 2 ) 



Therefore 



(10) 



Therefore, for translation invariant systems with Z^. sym- 
metry, the necessary conditions for a state be classical are 
either vanishing magnetization in the z spin direction or 
isotropy in the XY spin plane, i.e., 



or (atai) 



(11) 



Hence, the violation of the condition established in 
Eq. (jlip is sufficient for the nonclassicality of the quan- 
tum state. 



IV. FACTORIZED GROUND STATES IN THE 
XY MODEL 

As an example of a ^-symmetric system, let us con- 
sider a XY spin-1/2 chain, which is governed by the fol- 
lowing Hamiltonian 



H 



E 

i=0 



{I [(1 + 7K< + i + (1 - 7KVf +1 ] + ha\ 



(12) 

with N being the number of spins in the chain, cr™ the 
i-th spin Pauli operator in the direction m = x,y, z and 
periodic boundary conditions assumed. For the (dimen- 
sionless) anisotropy parameter 7 — ¥ 0, the model reduces 



to the XX model whereas for all the interval < 7 < 1 
it belongs to the Ising universality class, reducing to the 
transverse field Ising model at 7 — 1. The parameter h is 
associated with the external transverse magnetic field. In 
the thermodynamical limit, for A = J/h = 1, a second- 
order quantum critical line takes place for any < 7 < 1. 

The exact solution of the XY model is well known ~fl9l. 
[20j . The Hamiltonian (|12l) can be diagonalized via 
a Jordan- Wigner map followed by a Bogoliubov and 
Fourier transformation. Given the Zi symmetry, the two- 
spin reduced density operator can be written as 



P= 7 



I®I- 



3 

E 

i=l 



c 4 (/( 



(7° + cr" 



where, in terms of the parametrization given in Eq. ([7]), 
we have c\ = 2z + 2/, C2 = 2z — 2f, C3 = a + d — 2b, and 
C4 = a — d (with &x = 63 = 6). Due to the fact that the 
system is translationally invariant, the reduced state (fT3)) 
depends only on the distance between the spins. We will 
consider here nearest-neighbor pairs. The magnetization 
density G 2 = (a z ) and the two-point functions G^ +1 = 
( a a a ' t p~ 1 ) can b e directly obtained from the exact solution 
of the model [Hill. 

In the ordered ferromagnetic phase (A > 1), a finite 
magnetization (erf) (in the coupling direction) emerges, 
breaking the Z2 symmetry of the Hamiltonian for an infi- 
nite chain. This is possible because, at the critical point, 
the ground state becomes two-fold degenerated, exhibit- 
ing opposites values of (of) (denoted as \+m) and \—m)). 
Naturally, a superposition of the two ground states is 
also a ground state, which could then preserve the Z2 
symmetry. However, in the thermodynamic limit, these 
two degenerated ground states are not connected by local 
unitaries. Thus, to go from one to the other, one needs 
to flip all the spins at the same time, which is both a 
highly non-local and energy costly operation. Therefore, 
a small external perturbation will always pick up one of 
the two states, and the system will be "frozen" at it. The 
use of a superposition that preserves the symmetry could 
be justified in the case of a finite system, where the sys- 
tem, if highly isolated, could exist in such state. Another 
possibility would be to consider an equal mixture of the 
two ground states, which could be reasonable in the case 
one prepares a system at some temperature T and cools it 
down to the limit of T — > 0. This is a ^-symmetric state 
called thermal ground state. Note that, since the thermo- 
dynamic properties depend only on the energy spectrum 
(and not on the eigenstates), all choices of ground states 
will be equivalent. However, classical and quantum cor- 
relations do depend on the eigenstates and, therefore, the 
symmetric and broken choices of ground states may ex- 
hibit different pattern of correlations plM^ . 

Remarkably, the ground state with SSB exhibits a fac- 
torization (product state) point at 



7 2 + 



(14) 
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This point identifies a change in the behavior of the cor- 
relation functions decay, which pass from monotonically 
to an oscillatory decay [20]. Besides, it has been realized 
that, at this point, the product ground state is two-fold 
degenerated even for finite systems [22T - I241 ] . 

In order to apply the witness ||W|| in the XY model, we 
will first consider the case of the thermal (^-symmetric) 
ground state, for which the reduced density matrix is 
given by Eq. (|T3)) . For this case, a plot of ||W|| as a 
function of A is provided in Fig. [1] for 7 = 0.6. Observe 
that this plot shows that the only possible classical state 
appears for A = 0. Indeed, this point corresponds to 
an infinite transverse magnetic field applied, which leads 
the system to a product state, with all spins individually 
pointing in the z direction. Concerning the factorization 
point given in Eq. (|14l) , it is absent for the thermal state, 
since the statistical mixture of non- orthogonal product 
states exhibits nonvanishing quantum correlations [23j . 



Witness with Z^symmetry 
■ " Witness with SSB 
•— • Exact discord with SSB 

0.2 



0.15 




FIG. 1. (Color online) Witness ||W|| and exact value of QD 
for a pair of spins in both symmetric and broken ground states 
of the XY spin-1/2 infinite chain as a function of A for 7 = 0.6. 

On the other hand, if SSB is taken into account, which 
must be the case in the thermodynamic limit, this picture 
dramatically changes. Indeed, in the broken case, the 
two-spin reduced density matrix is more complicated and 
reads 



Pab 



( a p p f 

p b z q 

p z b q 

\ f q q d 



(15) 



where the functions p(A, 7) and g(A, 7) are given by 



P(A,7) 
e?(A,7) 



(at 



(16) 
(17) 



The evaluation of these functions have been detailed 
discussed in Ref. |2lj] . where tight lower and upper 
bounds for p(X, 7) and q(X, 7) are obtained. In our plots, 
either of such bounds essentially yields the same curves, 
which led us to keep the results produced with the lower 
bound. Concerning the results for ||W||, we can compare 
the plots for the thermal ground state and the ground 



state with SSB also in Fig. Q] Observe that, in the 
case with SSB, the only point for which ||W|| = is 
A = 1.25. Indeed, this corresponds to a classical state, 
as can be confirmed by the exact computation of the 
symmetric QD. Remarkably, this classical state is associ- 
ated with one of two doubly-degenerated factorized (fully 
product) ground state, which appears as a consequence 
of the Z2 SSB. This characterization of factorized ground 
states through ||W|| =0 can be numerically confirmed in 
Fig. [2~1 where ||W|| is plotted for a number of different 
values of the parameter 7. These unique classical points 
for the XY model are then clearly revealed by the wit- 
ness evaluation. Another aspect of the witness observed 
from Figs. Q] and [2~] is that it exhibits nonanalyticity in 
its first derivative at the second-order quantum critical 
point A = 1. This phenomenon occurs both for the sym- 
metric and broken ground states, which promotes ||W|| 
to a simple useful tool also to detect quantum phase tran- 
sitions. 



£ 0.15 




FIG. 2. (Color online) Witness ||W|| for a pair of spins in the 
ground state with SSB of the XY spin-1/2 infinite chain as a 
function of A for several values of the parameter 7. 



V. QUANTUM CORRELATED STATES IN THE 
ASHKIN-TELLER MODEL 



The Ashkin- Teller model has been introduced as a gen- 
eralization of the Ising spin-1/2 model to investigate the 
statitiscs of lattices with four-state interacting sites [25[ . 
It exhibits a rich phase diagram [2(| and has recently at- 
tracted a great deal of attention due to several applica- 
tions, e.g., nonabelian anyon models [27| and orbital cur- 
rent loops in Cu02-plaquettes of high-Tc cuprates [28$ . 
The Hamiltonian for the quantum Ashkin- Teller model 
in one-dimension for a chain with M sites is given by 



H AT = -J 



M 

E 

3=1 



t j t j 



M 

i=i 



^Z^Z 1 A —Z 2 

i 7+1 7 7 



-1), (18) 



5 



where er" and r" (a — x,y,z) are independent Pauli 
spin-1/2 operators, J is the exchange coupling constant, 
A and j3 are (dimensionless) parameters, and periodic 
boundary conditions (PBC) are adopted, i.e., <jf 1+1 = erf 
and Tj^ .J = t" (a = x,y,z). The Ashkin- Teller model 
is (g> .^-symmetric, with the Hamiltonian commuting 
with the parity operators 



and 



r 2 



(19) 



Therefore, the eigenspace of Hat can be decomposed into 
four disjoint sectors labelled by the eigenvalues of V\ and 
7^2 ■ We then numerically diagonalize Hat (via power 
method) to evaluate the witness ||W||, which provides a 
sufficient condition for the absence of Zi ® ^-symmetric 
classical states for finite values of (3 and A. We illus- 
trate this result by considering different configurations of 
multipartite states. More specifically, we choose a spin 
quartet, namely, a block of L = 4 particles described by 
spin operators ^ J ^Lj frkTki and respective generalization 
for a spin octet. For such cases, the absence of classical 
states for any finite A is displayed in Fig. [3] for (3 = 1 in 
a chain with N — 16 particles. This is confirmed in the 
inset of Fig. [3J where the exact behavior of global QD is 
exhibited. 
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FIG. 3. (Color online) Witness ||W|| as a function of A for 
several L-spin multipartite states in the Ashkin- Teller model 
for P = 1 in a chain with N = 16 spin particles. 

Moreover, besides sufficent for identifying no classical 
states in the thermal (Z 2 <S> ^-symmetric) ground state of 
the Ashkin- Teller chain, the witness ||W|| can also be ap- 
plied to identify second-order quantum phase transitions 
driven by the parameter (3, We illustrate this result by 
plotting || W|| for a quartet in Fig. [4] as a function of (3 for 
A = 3. From this plot, a second-order critical point can 
be readout through a minimum in the derivative of ||W| 
(see inset in Fig. [4}, which gets pronounced as the length 
of the chain is increased (we observe that the maximum 
appearing in the plot for d||W||/d/3 does not indicate any 
transition, since it does not exhibit pronunciation as we 



increase the size of the chain). This pronounced mini- 
mum, which is a typical precursor of the nonanalyticity 
at the thermodynamic limit, occurs for j3 — 0.61, which is 
compatible with the quantum phase diagram in Ref . [26j j . 
A further transition would be expected to occur near 
j3 = 1.5. This critical point is not apparent from ||W||, 
at least for chains up to N = 20 spins. Indeed, this sec- 
ond critical point is usually harder to characterize (see, 
e.g., characterization with entanglement in Ref. [2{| and 
with global QD in the upper insets of Fig. HI where pro- 
nunciation is subtle). 
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FIG. 4. (Color online) Witness ||W|| as a function of /3 for a 
quartet in the Ashkin- Teller model for A = 3 in chains with 
N = 8 and N = 20 spin particles. 



VI. CONCLUSION 

In summary, we have proposed a multipartite witness 
for nonclassical states based on their disturbance under 
local measurements. This witness provides a sufficient 
condition for nonclassicality that is in agreement with a 
nonvanishing global QD. As an illustration, we derived 
the necessary conditions for classicality in terms of corre- 
lation functions for ^-symmetric density operators and 
apply the witness for both the XY and Ashkin- Teller 
chains. Remarkably, the witness detects the factorized 
ground state of the XY chain as its unique classical state 
(for finite value of A and 7) as the SSB was taken into 
account. Then, in this example, ||W|| turns out to be 
not only sufficient but also necessary for nonclassicality. 
Moreover, the witness was also shown to be able to iden- 
tify second-order quantum phase transitions through a 
nonanalyticity in the first derivative of ||W|| for both XY 
and Ashkin- Teller chains. Even though classicality is rare 
in many-body microscopic systems |15| . it is rather use- 
ful, as a resource to both quantum information and con- 
densed matter applications, to characterize the existence 
of nonclassical correlation through simple and direct pro- 
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cedures. In particular, the typicality of states with van- 
ishing genuine n-partite correlations would be a possible 
application of our results. The witness ||W|| may consti- 
tute a helpful tool for such kind of characterization. 
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